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ABSTRACT 

We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 /im of the primary transit of the 
exoplanet GJ436b, by using blind source separation techniques. These observations are impor¬ 
tant to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed 
strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts 
on the possibility to extract conclusively an atmospheric signal; those analyses also reported dis- 
crepant resu lts, hence the necessity of this reanalysis. The method we used has been proposed in 
More llo ct alj ( 2014 1 to analyze 3.6 /iin transit light-curves of the hot Jupiter HD189733b; it per- 
formes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series 
read by individual pixels, from the same photometric observation. Our method only assumes the 
independence of instrumental and astrophysical signals, and therefore guarantees a higher degree 
of objectivity compared to parametric detrending techniques published in the literature. The 
datasets we analyzed in this paper represent a more challenging test compared to the previous 
ones. 

Contrary to previous results reported in the literature, our results (1) do not support any de¬ 
tectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable 
at the level ~10 -4 in the IR, and (3) the transit depth measurements at the two wavelengths are 
consistent within ler. We also (4) detect a possible transit duration variation (TDV) of ~80 s (2 
a significance level), that has not been pointed out in the literature, and (5) confirm no transit 
timing variations (TTVs) >30 s. 

Subject headings: methods: data analysis - techniques: photometric - planets and satellites: atmospheres 
- planets and satellites: individual(GJ436b) 


1. Introduction 


large majority of transiting exoplanets are “hot 
Jupiters”, i.e. planets with size similar to Jupiter 
orbiting very closely to their host star (semimajor 
axis ~ 0.01 — 0.5AU). Their typical surface tem- 


Transit spectroscopy and differential photome¬ 
try are largely used to investigate the composition 
and structure of exoplanetary atmospheres. The 


1 




peratures are > 1000K. 

GJ436b is a Neptune-sized planet orbiting 
around an M dwarf with radius ~0.46f?Q at 
a distance ~0.03 AU. This planet is interest¬ 
ing for several reasons. It is one of the small¬ 
est (radius ~4.3i?®) and coolest (~ 700K) ex- 


been measured (Gillon et al. 2007; Demory et al. 

2007: Alonso et al. 

2008s Coughlin et al. 2008s 

Caceres et al.l 2009: 

Doming et al. 2009; Pont et al. 

2009; Ballard et al 

20 ld 

Stevenson et al. 2010s 

Beaulieu et al. 201 

; Knutson et al. 2011. 20141). 


The primary transit depth is ~0.7%. Another pe¬ 
culiarity of GJ436b is its high orbital eccentricity 
(e ^0.1 6), inferred from ra dial velocity measure¬ 


ments (IManess et al.1 12 0071) and from secondary 


eclipse phasing (iPeming et al.l 1200911 . Both the 


physical and dynamical properties of GJ436b are 
debated in the literature. 

IManess et al. ( 2007 1 and Demorv et al. ( 2007 1 


investigated the origin of the high orbital ec¬ 
centricity of GJ436b, concluding that the cir¬ 
cularization timescale (^ 10 8 yr) is significantly 
smaller than the a ge of the system (> 6 xl 0 9 ). 


Maness et al.1 0200711 also found a long trend in 


radial velocity measurements; they suggested the 
presence of an external perturber on a wider orbit 
to explain both the high eccentricity of GJ436b 
and the lon g trend in radial velocity measure¬ 


ments. IRibas et al.l (120081 1 hypothesized a Super- 


Earth on a close orbit to explain those evidences, 
later retracte d. Transit tim i ng va ri ations (TTYs 


report ed by Alonso et al. ( 2008h : Caceres et al 


(l2009ll do n ot support any evidenc e of external 
perturbers. Stevenson et al.1 ( 2012ll claimed the 
possible detection of two nearby sub-Earth-sized 
exoplanets transiting in GJ436 system; accord¬ 
ing to the authors, the dynamic of the proposed 
system is consistent with the current non-TTV- 
detections. 

Based on multiwavelength i nfrare d eclipse mea¬ 
surements, Stevenson et al. ( 201Clh proposed a 
high CO-to-CH 4 ratio compared to thermochem¬ 
ical equilibrium models for hydrogen-dominated 
atmospheres. Their atmospheric model includes 
disequilibrium processes, such as vertical mixing 
and polymerization of met han to explain t he ob- 
served deficiency of CH 4 . iBeaulieu et al.l (1201111 
suggested strong CH 4 absorption at 3.6, 4.5, and 
8.0 /im Spitzer/IRAC passbands from primary 


transit observations, and their reanalysis of sec- 
ondar y eclipse data i s cons istent with this detec¬ 
tion. Knuts on et al.l |2011) measured significant 
time variations of the transit depths at the same 
wavelengths, which strongly affect the inferred 
transmission spectrum. They attributed such 
variations to the stellar activity and found that 
different results are obtainable depending on the 
observations considered. By rejecting those obser¬ 
vations that they believe to be most strongly af¬ 
fected by stellar activity, their final results support 
CO as the dominant carbon molecule, with very 
little, if any, CH 4 . More recent Hubble/WFC3 
observations in t he 1 . 2 — 1.6 urn wavel ength inter¬ 
val, analyzed by iKnutson et al.l ( 2014h . indicate a 
featureless transmission spectrum, which is con¬ 
sistent with relatively hydrogen-poor atmosphere 
with a high cloud or haze layer. 

In this paper we reanalyze four transit light- 
curves obtained with Spitzer/IRAC at 3.6 and 
4.5 /rm passbands (channels 1 and 2 of IRAC). 
We adopt a non-parametric data detrending tech¬ 
nique, based on Independent Component Analysis 
(ICA) applied to single pixel-light-curves, to en¬ 
sure a higher degree of objectivity. This method 
has proven to give robust results, when applied 
to the transits of the hot-Jupiter HD189733b 


obser ved with IRAC at 3.6 /im (IMorello et al 
Ml. We further test here the performance of 


this detrending technique with the more challeng¬ 
ing datasets of the Neptune-sized planet GJ436b, 
for which the transit depth is comparable with the 
amplitude of the instrumental pixel-phase signal, 
and the transit duration is very similar to the pe¬ 
riod of that signal. Additionally, we discuss the 
stellar and orbital stability of the GJ436 system, 
the repeatability of transit measurements, poten¬ 
tially affected by stellar, planet, and instrument 
variability, and the atmospheric contribution. We 
discuss the reliability of our results in light of other 
obse rvations reported in t h e literature, in partic¬ 
ular Beaulie u et al.1 ( 2011 1: Knutson et al.1 ( 2011 


20141) . 


2. Data Analysis 
2.1. Observations 

We analyze four photometric observations of 
GJ436b, which are part of the Spitzer program 
ID 50051. They include two 3.6 and two 4.5 /im 
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Table 1: Spitzer observations of primary transits of GJ436b. 


Obs. Number 

Detector 

Wavelength (pm) 

UT Date 

Orbit Number 

la 

IRAC, chi 

3.6 

2009 Jan 9 

234 

lb 

IRAC, chi 

3.6 

2009 Jan 28 

241 

2 a 

IRAC, cli2 

4.5 

2009 Jan 17 

237 

2 b 

IRAC, ch2 

4.5 

2009 Jan 31 

242 


primary transits as detailed in Tab. [T] Each ob¬ 
servation consists of 1829 exposures using IRAC’s 
sub-array mode, taken over 4.3 hr: 0.8 hr on the 
primary transit of the planet, the remaining 3.5 
hr before and after transit. The interval between 
consecutive exposures is 8.4 s. Each exposure in¬ 
cludes 64 consecutive frames integrated over 0.1 
s. We replaced the single frames of each exposure 
with their averages to reduce the random scatter, 
and the computational timt0. During an observa¬ 
tion, the centroid of the star GJ436 is stable to 
within one pixel. 


2.2. Detrending method, light-curve fit¬ 
ting and error bars 


Here we outline the main steps of the analysis, 
i.e. data detrending, light-curve fitting and esti¬ 
mating par ameter error bars. F urther details are 
reported in lMorello et al.l ( 20141 '). 


To detrend the transit signals from single ob¬ 
servations, we performed an ICA decomposition 
over selected pixel-liglit-curves, i.e. time series 
from individual pixels. We considered 5x5 ar¬ 
rays of pixels with the stellar centroids at their 
centers. In this way, we obtain a set of maxi¬ 
mally independent components: one of them is 
the transit signal, others may be instrumental sys- 
tematics and/or astrophysical signals. Observed 
light-curves are linear combinations of these inde¬ 
pendent components, the coefficients of the linear 
combinations can be calculated by fitting the out- 
of-transit parts. To estimate the transit signal in 
a robust way, the fit is performed on the out-of¬ 
transit of the relevant integral light-curve, i.e. the 
sum of the pixel-light-curves from the array used, 


1 Computational time is dominated by the time for transit 
fitting, which strongly depends on the algorithms (and set¬ 
tings) used. In our case, transit fitting time was of several 
hours, and it scales as 0(dN ), being d the number of free 
transit parameters, and N the data points. 


including all the non-transit components plus a 
constant term. The detrended transit signal is ob¬ 
tained by subtracting all the non-transit compo¬ 
nents, properly scaled by their fitting coefficients, 
from the integral light-curve. It is renormalized 
by the mean value on the out-of-transit, so that 
the out-of-transit level is unity. 


After the extractions of the detrended and nor¬ 
malized tr ansit time series , we m odelled them by 
using the iMandel fc Agoll ( 2002h analytical for¬ 
mulae. We originally assumed the orbital period, 
P, and t he epoch of the first transit, E tr , re¬ 
ported by ICaceres et al. ( 20091 ): the eccentricity, 
and the arg um ent of periastron, w, reported 


by Maness et al. ( 20071) : they are co nsistent with 


2004]; Gillon et al. 20071; Deming et al. 

200 $, and 

more accurate. We tested two differ 
quadratic limb darkening coefficients f 

ent sets of 
or the star 
by an At- 
a Phoenix 

(Howarth 2011b), 7 i and 73 , derived 

las (Kurucz 1970; Howarth 2011) and 

(Allard & Hauschildt 19951: Allard et al. 200lh 


models (see Sec. 12.3.11) . With these settings, we 
estimated the planet-to-star radii ratio, p = 
the orbital semimajor axis in units of stellar radii, 
ao = -g-, and inclination, i. First estimates were 
obtained through a Nelder-M ead optimization al¬ 


gorithm ( Laearias et al. 19981) : they were used as 


optimal starting points for an Adapt ive Metropo¬ 


lis a lgorithm with delayed rejection ( Haario et al 
2006h . generating chains of 20,000 values. Up¬ 


dated best estimates and (partial) error bars of the 
parameters, a par ,o , are the means and standard 
deviations of the relevant (gaussian distributed) 
sampled chains, respectively. The final parameter 
error bars are: 


®par 


&par, 0 


/ + rtcA 

V -0 


(1) 


ctq is the sampled likelihood variance, approxi¬ 
mately equal to the variance of the residuals for 
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the best transit model; cr\ CA is a term estimating 
the uncertainty associated to the ICA extraction 
(see App. IA..'tl for further details). 

For completeness, and for comparison with the 
literature, we also calculated the transit depth, p 2 , 
the impact pa rameter, b 1 and_the transit duration, 
T, where (see lFord et al. < 200Sh ) : 


b = Oq cos i 


1-e 2 


1 + e sin u> 


T = 


PVT^W V / T 3 e 2 

7rao 1 4- e sin to 


For a more thorough analysis, we performed 
other fits with different choices of the free param¬ 
eters, introducing a phase-shift parameter to con¬ 
sider possible timing error/variations, and simul¬ 
taneous fits on more than one multiple light-curves 
with some common free parameters. 

2.3. Application to observations 

Fig. [I] reports the raw “integral light-curves” 
observed. The main systematic effect for IRAC 
channels 1 and 2 observations is an almost regular 
undulation with period ^3000 s, so-called pixel- 
phase effect, because it depends on the relative po¬ 
sition of the_source centro id with respect to a p ixel 
cente r ( Fazio et, al. l20f)4 iMorales- Calderon et al. 
20061 ). This effect is particularly difficult to de¬ 


trend from these datasets because its timescale is 
similar to the transit duration, and its amplitude 
is comparable to the transit depth. Recently, a 
time depen dence of the pixel-phas e effect has been 


sugge sted (Stevenson et al. 

201 i f. 


201ft Beaulieu et al 


If considering the whole datasets for detrend¬ 
ing, our ICA algorithm is able to remove most of 
the non-flatness on the out-of-transits, and visi¬ 
bly improve the in-transit shapes, but some vis¬ 
ible issues remain (see Fig. usd. We noted that 
results improve significantly if rejecting a number 
of data points from the beginning of each obser¬ 
vation. It is discussed on a statistical basis in 
Sec. 13.11 A possible explanation is that first data 
points contain a long-tail variatio n until stabiliza¬ 


tion of the instruments ( Fazio et al. (20041), 


see 


also App. EH) ; this is not a crucial point for the 
data analysis. In the rest of this paper we dis¬ 
cuss the results obtained after rejecting the first 
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mary transit observations. Data points on the left 
of black vertical lines have been discarded for the 
analysis. Note that the transit depth is compara¬ 
ble with the amplitude of systematics. 
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Fig. 2.— Central pixel-light-curves of the four pri¬ 
mary transit observations. Data points on the left 
of black vertical lines have been discarded for the 
analysis. Note that pixel-light-curves are domi¬ 
nated by systematics, and the transit signal is not 
visible by eye (but it is present, as proven by ICA 
retrieval). 

450 exposures from each observation, correspond¬ 
ing to ^3780 s, for which the ICA performances 
are optimal. It is worth to point out that different 
choices (including no data rejections) give consis¬ 
tent results (within 1 cr), with larger or similar 
error bars. 

2.3.1. Limb darkening coefficients 

Tab. [2]reports the quadratic limb darkening co¬ 
efficients used at 3.6 and 4.5 pm IRAC passbands. 
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Both the Atlas and Phoenix model s are com p ute d 


with T e fj = 3680 K, lo gq = 4 .78 (T orres 


and solar abundances (jAsnlund et, al.ll2009Tl . 


2009), 


Table 2: Quadratic limb darkening coefficients 
computed by Atlas and Phoenix stellar models for 
IRAC 3.6 and 4.5 pm bands. 


Atlas 

3.6 /im 

4.5 pm 

7i 

5.489 x 10 -2 

1.331 x 10 -2 

72 

3.0653 x 10" 1 

2.8396 x 10 1 

Phoenix 

3.6 /im 

4.5 pm 

7i 

3.87 x 10- 3 

3.27 x 10 -3 

72 

2.3615 x lO' 1 

1.8193 x 10 1 


3. Results 

3.1. Tests of pixel-phase correlations 

To investigate the effectiveness of the data de¬ 
trending we measure the correlations of the sig¬ 
nals with the pixel-phase position, before and after 
the corrections. We refer to the Pearson product- 
moment correlation coefficient (PCC), defined as: 

PCC = mv{X - Y) (4) 

ax cry 

where cov(X, Y) is the covariance of the signals X 
and y, ct x and ay are the standard deviations. 
In this context, X and Y are temporal series of 
fluxes and pixel-phases. The PCCs are measured 
over three intervals, i.e. pre-, in-, and post-transit, 
where the astrophysical signals are expected to be 
almost flatH. In general -1<PCC<+1, where +1 
is total positive correlation, -1 is total negative 
correlation, and 0 is no correlation. Fig. [3]reports 
the temporal series of pixel-phases. Fig. [4] reports 
the values of the PCCs measured on the pre-, in-, 
and post-transit for each observation, uncorrected, 
corrected without and with pre-transit truncation. 
The original data are strongly anticorrelated with 
the pixel-phase, with PCC<-0.9 for channel 1, and 


2 We used the following definitions: pre-transit (fp < - 

0.0082); in-transit (-0.00433 < (j> <0.00416); post-transit 
( (f> >0.0082). These have been decided so that all the tran¬ 
sit models obtained during the analysis, modified with no 
limb darkening, are flat in these three intervals. We checked 

that other reasonable choices of the limits do not affect this 
analysis. 


PCC~-0.7 for channel 2. After the ICA detrend¬ 
ing including all the data, these correlations are 
significatively reduced (|PCC| <0.3). If we remove 
the first 450 data points, the ICA detrending gen¬ 
erally performs significantly better (|PCC| ~10 -3 - 
7xl0 -2 ). Fig. [5] reports the level of significance 
of the residual correlations in the detrended data, 
calculated with a permutation test. When we re¬ 
ject the first 450 data points, the residual correla¬ 
tions in the detrended data are below 1.5 er, except 
for Obs. 2a, for which the residual correlation is 
higher in any case. The residual correlations with¬ 
out the cut of the first 450 data points are larger, 
i.e. >4 cr in the post transit, with the exception 
of Obs. lb, for which the residual correlations are 
below 1 a in any case. 




Fig. 3.— Time series of the pixel-phase values 
for the four observations. Data points on the 
left of black vertical lines have been discarded 
for the analysis; dashed green lines delimit the 
ends of pre-transits and the begins of post-transits; 
dashed red lines delimit the in-transits. 


3.2. Fitting p, ao and i 

Fig. [ 6 ] reports the detrended light-curves, 
binned over 7 points, with the relative best transit 
models, and the residuals. The transit models in 
Fig. [ 6 ] are computed with 71 and 72 Phoenix coef¬ 
ficients. Analogous transit models computed with 
71 and 72 Atlas coefficients are very similar, with 
average standard deviations <1.9xl0 -5 , and max¬ 
imum discrepacies <K) _4 . Discrepancies between 
the transit models and the detrended light-curves 
are at the level ~2.0xl0 -4 for IRAC channel 1, 
and ~2.6—2.9x 10 ~ 4 for IRAC channel 2, there- 
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3.2.1. Combining observations 
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Fig. 4.— PCCs between fluxes and pixel-phases 
for pre-, in-, and post-transits of the four light- 
curves; (blu circles) raw data, (green rightwards 
triangles) ICA detrended data with no rejections, 
(red upwards triangles) ICA detrended data after 
rejecting the first 450 points. 



Fig. 5.— Significance level of correlation between 
fluxes and pixel-phases for the four observations; 
(green rightwards triangles) ICA detrended data 
with no rejections, (red upwards triangles) ICA 
detrended data after rejecting the first 450 points. 


fore it is not possible to distinguish between Atlas 
and Phoenix models from the data. Best param¬ 
eter results and error bars are reported in Fig. 0 
in Tab. H] and [5] Atlas and Phoenix stellar mod¬ 
els lead to two systematically different parameter 
sets, but within the error bars. All the parameters 
from different observations are comparable within 
la, even neglecting the detrending errors (ct/ca); 
except the transit durations for Obs la and lb. 
This is discussed in the following sections. 


We performed two couples of simultaneos fits, 
one for the 3.6 (im and one for the 4.5 pm light- 
curves, with Atlas and Phoenix limb darkening 
coefficients, assuming common orbital parameters 
(ao and i), and potentially different transit depths 
(p), in order to cancel the effects of parameter in¬ 
tercorrelations. The assumption that orbital pa¬ 
rameters are the same during each observation 
is very reliable, because they are sparse over a 
short period of time (less than 1 month, 9 plan¬ 
etary orbital periods), so that variations due to 
relativistic effects, external perturbers or tidal ef¬ 
fects, would be very sma ll compared to the error 
bars_ (I Alonsoi et alJ 12008 : Ijordkn fe Bakosl 12008 : 
Pal fe Kocsis 20081) . 


The results of these combined fits are reported 
in Fig. E El in Tab. [HI and 0 The 4.5 pm tran¬ 
sit depths become identical, with an intermedi¬ 
ate value between the two determined with sep¬ 
arate fits; the 3.6 pm transit depths slightly di¬ 
verges, but their separation is still less than ler. 
The standard deviations of residuals between the 
detrended light-curves and the transit models in¬ 
crease of ~2—3xl0~ 6 for Obs 2a and 2b (neg¬ 
ligible), and of ~7—8xl0 -6 for Obs la and lb 
(comparable with the op uncertainties). The as¬ 
sumption of common orbital parameters for Obs 
2a and 2b may be valid, being the consequent 
transit models as good as the individually fitted 
ones. Being the transit depths also identical, the 
two light-curves are very well approximated by the 
same transit model. The original discrepancies 
between the two sets of transit parameters were 
enlarged by their intercorrelations. The same as¬ 
sumption for Obs la and lb lead to worse transit 
models and more divergent transit depths, but, in 
both cases, not dramatically. 


3.3. Timing variations 

We performed transit model-fits with a free 
phase-shift in addition to p , ao, and i, in order 
to investigate the effect of possible timing varia¬ 
tions. Fig. [lU]reports the time-shifts obtained. No 
evidence of timing variation have been detected, 
with upper limits <30 s. Both Atlas and Phoenix 
stellar models lead to the same shifts. Other pa¬ 
rameter estimates are not affected. 
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Fig. 6.- Left panels: (blue) detrended light- 
curves for the four observations with (red) best 
transit models overplotted, binned over 7 points; 
best transit models are calculated with p, a o, and 
i as free parameters, and Phoenix quadratic limb 
darkening coefficients (see Sec. 12.3.11) . Right pan¬ 
els: Residuals between detrended light-curves and 
best transit models; black horizontal dashed lines 
indicate the standard deviations of residuals. 
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Fig. 7. From top to bottom: Comparisons of 
the parameters p , a 0 , and i (left side), p 2 , b , and 
T (right side), obtained in this paper with Atlas 
stellar model (cyan, empty c ircles), Phoenix stellar 
model (blue, full circle s), in Knutson et al.1 (201 ll) 
(red triangles), and in Beaulieu et al.1 ( 201lh (yel¬ 
low squares). 
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Fig. 8.— Comparisons of the parameters p and p 2 , 
obtained in this paper, with common orbital pa¬ 
rameters for observations at the same wavelength, 
and Atlas or Phoenix stellar models. 
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Fig. 9.— Comparisons of the parameters b and T, 
obtained in this paper, with common orbital pa¬ 
rameters for observations at the same wavelength, 
and Atlas or Phoenix stellar models. 


Fig. 10.- Best fitted time-shifts of the mid¬ 
transit with respect to the periodically predicted 
mid-transit times, assuming Atlas or Phoenix 
stellar models. Note that results are model- 
independent, because mid-transit time will not 
correlate with limb darkening parameters, or any 
other physical parameters. 



















































































4. Discussion 

4.1. Comparing observations 

Fig. im and IT21 report the superpositions of 3.6 
and 4.5 //in liglit-curves respectively, and the resid¬ 
uals. In both cases the mean value of the in-transit 
residuals is small (;$5xl0 -5 ), but the transit lb 
is clearly longer than transit la, as measured by 
transit duration (T) parameters. As T is function 
of the orbital parameters and stellar model, this 
is the reason why simultaneous fits with common 
orbital parameters and stellar model do not be¬ 
have very well. We also note that the ingresses of 
transits 2a and 2b have different slopes. 
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Fig. 11. Top panel: detrended light-curves for 
Obs la (blue), and for Obs lb (green). Bottom 
panel: Residuals between the two observations. 
Black dotted lines delimit the out-of-transit, red 
dotted lines delimit the in-transit (as defined in 
Sec. ED - 
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13.2.11 Simultaneous fits over the four observations 
with common orbital parameters do not add any 
information. 


4.2. Comparison with previous analyses of 
the same observations 


Fig. [3 reports the parameter values obtained 
in this paper for the indivi dual observations wit h 
the analogues r eport ed by iBeaulieu et al.1 (1201111 ; 


Knutson et al.l ((2011). Our results suggest a con¬ 
stant value of the transit depth (largely within la) 
both between the 3.6 and 4.5 um observations, and 
for the two wavelengths. iKnutson et al .l ( 201l( ) re¬ 
port variations of the transit depth with a 3.4tr 
significance between the two epochs at 3.6 /zm, 
and 2.1cr at 4.5 /zm, which they att ributed to stel¬ 
lar activity. IBeaulieu et al. ( 201 1) also obtained 
significant differences between different epochs at 
the same wavelength, but they attributed such 
discrepancies to an unfavorable transit-systematic 
phasing, then they discarded those epochs from 
the analysis. Our error bars are generally compa¬ 
rable to the ones reported in both previous papers, 
but in some cases they are larger up to a factor ~2. 
This is not surprising, because we are not making 
any prior assumptions about the si gnals, to guar¬ 
antee a high degree of objectivity (Morello et al.l 
2014 : Waldmann 20121 : IWaldmann et al. 2013h . 
We conclude that our detrending method lead to 
more robust results than the previous ones in the 
literature, and they show no evidence of stellar ac¬ 
tivity variations at ~10~ 4 photometric level. Re¬ 
cent results f rom Hubble /WFC3 observations at 


1.2—1.6 //m (Knutson et al.112014 ) also show no 


significant transit depth variations over four ob¬ 
servations in about 2 months. 


Fig. 12.- Top panel: detrended light-curves for 
Obs 2a (blue), and for Obs 2b (green). Bottom 
panel: Residuals between the two observations. 
Black dotted lines delimit the out-of-transit, red 
dotted lines delimit the in-transit (as defined in 
Sec. 13.11) . 

The difference between p 2 values at the two 
wavelengths is (on average) ~5xl0 -5 , then there 
is no evidence of differences in planetary atmo¬ 
sphere’s absorption at the two wavelengths. 

The orbital parameters at the two wavelengths 
are also comparable, as detailed in Sec. 13.21 and 


4.3. Comparison with other observations 


Fig. [13] compares our estimated transit depth 
values at 3.6 and 4.5 //m (averaged over the obser¬ 
vations at the same wavelen gth) with the most re - 
cent results at 1.2—1.6 fi m ([Knutson et al.l|20 14). 
The resulting spectrum is featureless within the er¬ 
ror bars. However, when comparing transit depth 
measurements at different wavelengths, we should 
ensure that the GJ436 system is uniformly mod¬ 
elled, i.e. same stellar model and orbital parame¬ 
ters. A uniform multiwavelength reanalysis is re¬ 
quired to confirm this result, and investigate po¬ 
tential small features. 
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Fig. 13.— Transit depth values obtained in this 
paper at 3.6 and 4.5 u m (blue circles): valu es at 
1.2—1.6 pm reported bv lKnutson et al. ( 2014h (red 
triangles). 


Our non-detection of TTVs higher than ~30 s 


in the infrared 

Alonso et al. 20081 Caceres et al. 

2009; 

Pont et al. 

2009; 

Ballard et al. 2010t Knutson et al. 

2011, 

20141. We measured a significant TDV (W80 


s) between Obs la and lb; we did not find any 
study of TDVs for GJ436b in the literatu r e, but 


injecting parameters from [Knutson et al. (2011) 
into our Eq. [3] we obtain a similar trend for 
the same observations. More observations are 
required to investigate the cause of the apparent 
TDV between Obs la and lb, whether it is due to 
a perturber (as currently required to explain the 
high orbital eccentricity), a stellar phenomenon, 
or something else. 


nificance level), but no significant TTVs; more 
measurements are required to investigate the pos¬ 
sible presence of a perturber, and its nature. Also, 
more uniform analyses at other wavelengths are 
required to get a more reliable transmission spec¬ 
trum. 

G. Morello is funded by UCL Perren/Impact 
scholarship (CJ4M/CJ0T). I. P. Waldmann is 
funded by the European Research Council Grant 
“Exolights”. G. Tinetti is funded by the Royal 
Society. G. Micela is supported by “Progetto Pre- 
miale - A Way to Other Worlds”, funded by Italian 
Minister for University and Scientific Research”. 


5. Conclusions 


We have applied a blind signa l-source se pa¬ 


ration method, firstly proposed by iMorello et al 


(12014 ). to analyze other photometric data of pri¬ 
mary transits of an exoplanet, and extending its 
validity to the Spitzer/IRAC 4.5 (im band. These 
datasets were more challenging to analyze, be¬ 
cause of the lower transit depth, comparable with 
the amplitude of the instrumental pixel-phase sig¬ 
nal, the transit duration, very similar to the period 
of said signal, and possible stellar variability. 

We obtain consistent results between transits 
at different epochs, ruling out stellar activity vari¬ 
ations within ~10~ 4 photometric level. We do 
not detect any significant difference for the transit 
depth at 3.6 and 4.5 /mi, neither with the re¬ 


cent measurements at 1.2—1.6 /mi (Knutson et al 


20141) . supporting the hypothesis of a flat trans¬ 
mission spectrum. We measure a TDV of 80 s 
between transits separated by 7 orbits (2 er sig- 
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A. ICA 


A.l. Rationale 

ICA is a special case of “blind source separation” technique, i.e. it aims to separate original source signals 
from observations with minimal assumptions. The assumptions for standard ICA [f] are: 


1. the source signals are statistically independent; 

2. observations are linear mixtures of the source signals; 

3. the number of distinguishable observations is not smaller than the number of sources. 


The first condition is easily verified if the signals have different origins, i.e. the astronomical target, 
other background objects, and the instruments. Addit ionally, some studies found that ICA algorithm s can 
separate also signals that are not exactly independent ( Hvvarinen fc Oia 20001 : Hvvarinen et al.l|2001 b 
The second condition is more questionable, given that some instrumental systematics might be multi¬ 
plicative rather than additive. Alternative ICA algorithms consider non-linear mixing of the source signals, 
but some additional in formation are r equired to perform the separation, and, in general, there is not a 
unique solution (Hvvarinen et alj [2001l) . Based on the following evidences, we found that, for Spitzer/IRAC 
light-curves, the classic assumption of linear mixing leads to reliable and robust results: 


detrended light-curves pre sent a low level of residual scatter, compared to the literature (IBeaulieu et al 
l201ll iKnutson et al. 2011 ): 


• planetary and stellar parameters measured at different epochs are consistent (this is not a necessary 
condition); 

• non-transit components have the same characteristics, e.g. periodicity and amplitude, of known in¬ 
strumental systematics. 


This seems to be i n contrast with the standard (empirical) pixel-phase effect method used to detrend Spitzer 


data ( Fazio et al. 20041 ): flux measurements are correlated with the position of the centroid on a pixel, the 


cause of this is assumed to be an intra-pixel sensitivity variation, hence the systematics model is multiplied 
to the astrophysical signal. We are now investigating this question through simulated observations (Morello 
et al., in prep.); we report here our preliminary results: 


• either inter- and intra-pixel effects (or both) can originate systematics similar to the ones observed in 
Spitzer; 

• inter-pixel effects are additive, as in our ICA model; 

• intra-pixel effects are not additive, but the ICA algorithm is still able to significantly reduce their 
presence in the light-curves (our simulations currently indicate a reduction by a factor of 7 for the 
amplitudes of systematic components from an original 3.5xl0~ 3 photometric level, outperforming the 
pixel-phase method by a factor of 2.3-3.3). 


The third condition is case-dependent, since the number of components is not known a priori, and the 
number of pixels is limited by the width of the PSF. Also, if all the pixels contain the same systematic signals 
with the same weights relative to the astrophysical signal, the pixel-light-curves would not be distinguishable, 
and separation would be impossibile. Given the results obtained, we infer that we have a sufficient number 
of distinguishable pixel-light-curves to detrend our signals up to a 2xl0~ 4 photometric precision. 

3 There are some variants of ICA which do include prior information, e.g. Barriga ct al. (2011); Igual ct al. (2002); Stone et al. 

120021) . 
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A.2. Performances of MULTICOMBI algorithm 


In this section, we discuss the ability of ICA to se parate d i fferent kinds of signals. It depends on the 
particular algorithm used, in our case MULTICOMBI ( Tichavskv et al. 2008 ). MULTICO MBI is a powerful 
tool, that optimally mixes two complementa ry algorithms, i.e. EFICA ( Koldovskv et al.ll2006l) . designed to 
separate non-gaussian signals, and WASOBI ( Yeredor 2000lh specialized to separate gaussian auto-regressive 
and time-correlated components. 

We tested MULTICOMBI performances with simulated observations of planetary transits affected by a 
large variety of systematic signals, including non-stationary signals with changing frequencies and amplitudes, 
sudden change points, transient behaviours, and long-term monotonic drifts. In all cases the algorithm 
successfully dentrended the systematic components, except monotonic drifts (see Fig. [14] and [15] for some 
examples). Spitzer/IRAC lightcurves typically start with a drift before stabilization; this would explain the 
improved performance if rejecting part of the earlier data points. 

A more detailed analysis of the performances of the ICA detrending method adopted in this paper for 
different instrument cases, i.e. inter- or intra-pixel effects, amplitude and frequency of the systematics, 
temporal structures, non-stationarity, individual pixel peculiarities, is ongoing (Morcllo et ah, in prep.). 

The calculation time for a MULTICOMBI session is, in our cases less than 3 s. We are considering 
a number of signals d =25, and N~10 3 data points. For non-binned datasets, i.e. N ~10 5 data points, 
the calculation time varies in the range 45-105 s, then it is, in general, case dependent. We measured the 
computation time for different val ues o f d and N , a nd found that the algorithm complexity scales as 0(d 2 N ), 
as predicted for EFICA in Koldovskv et al. ( 2006 ). 
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Fig. 14.- Left panels: simulations of normalized raw ligh-curves with different systematic effects, due to 
pointing jitter and different pixel responses. Jitter time series are, from top to bottom: sinusoidal, periodic 
sawtooth (smoothed), successive smooth sawtooths with increasing frequency, and sudden shift. Right panels: 
Correspondent detrended light-curves with pixel-ICA method. Retrieved parameters are consistent with the 
values adopted to generate the raw light-curves. 
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Fig. 15.' Left panels: simulations of normalized raw ligh-curves with different systematic effects, due to 
pointing jitter and strong intra-pixel variations. Jitter time series are: (top) sinusoidal, (bottom) successive 
smooth sawtooths with increasing frequency. Right panels: Correspondent detrended light-curves with (blue) 
pixel-ICA method, and (red, dashed) traditional pixel-phase method. Note that ICA method outperforms 
the pixel-phase method: lower residual systematics and retrieved parameters closer to the true values. 
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A.3. ICA errors 


If ICA were able to separate the original source signals perfectly, the parameter error bars would be fully 
determined by the residual scatter on the detrended light-curve. In general, we expect this not to be the 
case, since any detrending method would introduce some bias in the parameter estim ates. We mode l su ch 
unknown bias as an additive uncertainty, eric a, hr the time series, leading to Eq. |T] iMorello et al. ( 2014 ) 
report the following formula for (Jica : 


a ICA — /“ ( C^ISRj + (Tntc-fit 


(Al) 


where ISR is the so-called Interference-to-Signal-Ratio matrix, Oj are the coefficients of the non-transit- 
components, m is their number, a n t c -fit is the standard deviation of residuals from the referent raw light- 
curve, out of the transit, / is the normalising factor for the detrended light-curve. The sum on the left takes 
into account the precision of the components extracted by the algorithm; a n t c -fit indicates how well the 
linear combination of components approximates the out-of-transit. Note that, while the first term increases 
with the number of components considered (see Sec. lT2l) . the second term decreases. The optimal strategy is 
to remove all the extracted non-transit components from the raw light-curve, though many resul ts obt ai ned 
by removing the most significant components (to be determined) are almost identical (see Morello et al 
(120141) . Sec. 2.5.2). 


MULTICOMBI code produces two Interference-to-Signal-Ratio m atrices , ISR/ 7 . asso ciated to the algo¬ 
rithm EFICA, and ISR 1174 , associated to the algorithm WASOBI. In Morello et al. ( 2014 ) we estimated the 
global ISR as the arithmetic mean of ISR^ and ISR 11 A ; this a very conservative estimate, which does 
not take into account the outperforming separation capabilities of MULTICOMBI compared to EFICA and 
WASOBI; here we suggest a more appropriate definition: 


ISRij = mm 


in (lSRff,ISR 


1,0 


M ) 


(A2) 


In the cases analyzed in iMorello et al. ( 2014 ) the contributions of the ISR terms were ~10% of the total 
error bars, then adopting the new definition of ISR would not modify the results significantly. Here, we 
find that the ISR contributions to the error bars are comparable with the other terms, probably because 
some instrumental systematics and the transit signals have similar timescales and amplitudes, making the 
separation more uncertain. Tab. [3] reports the values of <jica obtained for each observations, with ISR 
calculated according to Eq. IA21 and according to the arithmetic mean definition. 


Table 3: Estimated ct/ca values for the four observations (Eq. IAll and lX2l) . and worst case values according 
to a more conservative estimate of ISR (arithmetic mean of ISR^ and ISR" 4 ). 


Obs. number 

(tic a 

aicA (max) 

la 

4.24x 10 -4 

5.97xl0 -4 

lb 

3.41xl0” 4 

4.52xl0 -4 

2a 

2.79xl0 -4 

3.07x 10 -4 

2b 

3.97xl0~ 4 

6.42x 10 -4 


The error bars obtained in this p ape r with the definition in Eq. IA2l are consistent with results from the 


same tests of robustness reported in Morello et al. ( 20141) . i.e. different number of components, pixel arrays, 


and partial time series. For completeness, tables in App. [D] report the error bars obtained from residual 
scatter only, including uica with ISR defined as in Eq. IA21 and with the previous definition of ISR. The 
relative difference between error bars obtained with the two definitions of ISR is 0—30%, which, in general, 
may be important for atmospheric characterization, but it does not affect the conclusions obtained in this 
paper. 
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B. Full datasets analysis 

Fig. [16] reports the detrended light-curves obtained using the whole datasets. As stated in Sec. 12.31 
the systematics originally present in the raw light-curves are greatly reduced (see Fig. [T]), but there are 
still visible trends at all phases of the transits. Quantitative measurements of these trends, based on the 
correlations with the pixel-phase, have been discussed in Sec. 13.11 together with the trends obtained by 
rejecting the first 450 data points before processing. Also, standard deviation of residuals between the light- 
curves and the transit models are larger for the cases with no preliminary data rejection (but for Obs lb 
the difference is not significant). Removing the first 450 data points after ICA processing may reduce the 
scatter, but it does not improve the reliability of the light-curves, because trends are present at all times. 

Results extracted from these worse-quality light-curves are less robust, but they are consistent with 
accepted results within 1 a. Also, error bars are largely underestimated if neglecting the uncertainties due 
to the detrending process. It may induce to erroneously detect inter-epoch transit depth variations. 


C. Alternative transit model-fits 
C.l. Free limb darkening 

We performed transit model-fits with one free limb darkening parameter (linear or quadratic) in addition 
to the other free parameters (p, a o, and i). The standard deviations of residuals between the detrended 
light-curves and the transit models do not change, among the models obtained with Atlas, Phoenix, free 
linear and free quadratic limb darkening coefficients. The parameter error bars are larger by factors in the 
range 1—3 for the free quadratic case, even much larger for the free linear case. Best parameter estimates 
may be more affected by intercorrelations. The pure quadratic limb darkening is a better approximation 
of the real case, because in both Atlas and Phoenix models the quadratic coefficients are greater than the 
linear ones. Fig. [T3 reports the estimates for the quadratic limb darkening coefficients, 72 . 

Note that: 

1 . they are comparable (within lcr) with the theoretical values; 

2. error bars are larger than the differences between Atlas and Phoenix values; 

3. error bars do not allow to distinguish the values at 3.6 and 4.5 pm. 

Interestingly, the best estimate of 72 for Obs lb is the most distant from the other values (but within 
lcr), and the correspondent transit duration is now equal for Obs la and lb. Although the value 72 —0 is 
not reliable, it is important to note that the measured transit duration, as defined in Eq. [3] depends on the 
stellar intensity distribution and limb darkening model adopted, then: 

• stellar variability is a possible cause for observed TDVs; 

• TDVs measured from observations at different wavelengths must be taken carefully. 

C.2. Free eccentricity 

We performed transit model-fits with free eccentricity, e, in addition to the other free parameters (p, ao, 
and *). Eccentricity has not a great impact on the transit models: best estimates of the other parameters 
do not change significantly with respect to previous ones (with e = 0.16), although best estimates for the 
eccentricity for different observations varies over a large range (0.08—0.22). Also, residuals between light- 
curves and models are not affected. 
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Fig. 16.— Left panels: (blue) detrended light-curves for the four observations without data rejection, (red) 
best transit models overplotted, binned over 7 points; best transit models are calculated with p, ag, and i as 
free parameters, and Phoenix quadratic limb darkening coefficients (see Sec. 12. 3. 111 . Right panels: Residuals 
between detrended light-curves and best transit models; black horizontal dashed lines indicate the standard 
deviations of residuals. 

C.3. Free argument of periastron 

We performed transit model-fits with free argument of periastron, w, in addition to the other free param¬ 
eters (p, ag, and i). They do not constrain u> very well (<r w ~24—34); a o and i are strictly correlated with 
w, and their error bars are ~3 times larger than ones obtained with uj fixed. Also, the distributions of w, 
ao, and i, are asymmetric, because best w values are close to the edge of the range of admissible values. It 
is interesting to note that other parameters, such as p, 6, and T, are not affected by w degeneracies: their 
posterior distributions are indistinguishable from the ones obtained with oj fixed. 
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Atlas, IRAC 1 
Atlas, IRAC 2 
--Phoenix, IRAC 1 
--Phoenix, IRAC 2 


Fig. 17. Best fitted values of the quadratic limb darkening coefficient (72) for the four observations, 
assuming 71 = 0 . 

C.4. Free phase-shift 

We performed transit model-fits with a free phase-shift in addition to other free parameters, in order to 
investigate the effect of possible timing variations. Results are very similar for each fitting configuration, 
with no evidence of timing variation below ~30 s, as stated in Sec. 13.31 The free phase-shift does not affect 
other parameter estimates, except in the cases with free argument of periastron (w): timing variations and 
orbit precession are highly correlated. 

C.5. Fitting p , a 0 , and i with components’ coefficients 

We experimented an alternative method to estimate the transit parameters and the coefficients of the 
independent components simultaneously, by modeling the raw lightcurves as linear combinations of the 
components plus a transit model. In this way, we can investigate possible correlations between transit 
parameters and mixing coefficients, and test the stability of an ICA model over the whole observation. If 
results were significantly different than the ones obtained by estimating the mixing coefficients on the out-of- 
transit only, it would indicate that something different has happened during the transit, either astropliysical 
or instrumental in nature. For the datasets analyzed in this paper, results are consistent with accepted 
values within 1 a. It is worth to note that partial error bars, a part 0 , obtained from the MCMCs are similar 
to the ones obtained with p, ao, and i only free parameters (in some cases even smaller). This indicates that 
the main cause of uncertainty attributed to the detrending method is not given by the mixing coefficients, 
but the intrinsic errors on the components extracted. 



la 


1b 


2a 


2t 


D. Tables 
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Table 4: Transit parameter estimates for the four observations, by fitting p , a o, and i as free parameters, 
with ATLAS quadratic limb darkening coefficients (Tab. [2]). We report the partial error bars obtained by 
the residuals, the final error bars, and the worst case error bars (Eq. [lj Sec. 12.21 and App. IA.3D . 


Obs. number 

Parameters 

Best values 

l-fj errors 

(residual scatter only) 

1-a errors 
(ICA) 

1-a errors 
(ICA worst case) 

la 

P 

do 

i 

P 2 

b 

T( s) 

do 

0.0834 

14.24 

86.63 

0.00696 

0.836 

2835 

2.01x 10~ 4 

3x 10~ 4 

0.20 

0.07 

5x 10~ 5 

0.005 

10 

O.lOx 10~ 4 

7xl0 -4 

0.46 

0.16 

l.lxlO -4 

0.012 

24 

9x 10~ 4 

0.61 

0.21 

1.5x 10 -4 

0.016 

33 


P 

0.0836 

3xl0~ 4 

6xl0 -4 

7xl0“ 4 


do 

13.93 

0.19 

0.39 

0.48 


i 

86.57 

0.07 

0.14 

0.17 

lb 

p 2 

0.00699 

5xl0' 5 

9xl0 -5 

1.2xl0 -4 


b 

0.834 

0.005 

0.011 

0.013 


T(s) 

2916 

10 

20 

25 


do 

1.96x 10~ 4 

O.lOx 10 -4 




P 

0.0833 

3xl0' 4 

Cn 

X 

o 

1 

5x 10 -4 


do 

14.24 

0.25 

0.37 

0.39 


i 

86.66 

0.09 

0.13 

0.14 

2a 

p 2 

0.00694 

6xl0' 5 

00 

X 

o 

1 

or 

9x 10 -5 


b 

0.831 

0.007 

0.010 

0.011 


T(s) 

2879 

13 

19 

21 


do 

2.56x 10 -4 

0.13x 10 -4 




P 

0.0842 

4xl0' 4 

6xl0 -4 

9x 10“ 4 


do 

13.38 

0.23 

0.40 

0.57 


i 

86.34 

0.09 

0.15 

0.21 

2b 

p 2 

0.00709 

6x 10~ 5 

l.lxlO -4 

1.6x 10 -4 


b 

0.854 

0.006 

0.010 

0.014 


T( s) 

2860 

16 

27 

39 


do 

2.87x 10~ 4 

0.15x 10 -4 
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Table 5: Transit parameter estimates for the four observations, by fitting p , a o, and i as free parameters, 
with PHOENIX quadratic limb darkening coefficients (Tab. [2]). We report the partial error bars obtained 
by the residuals, the final error bars, and the worst case error bars (Eq. [TJ Sec. 12.21 and App. IA.3I) . 


Obs. number 

Parameters 

Best values 

l-fj errors 

(residual scatter only) 

1-a errors 
(ICA) 

1-cr errors 
(ICA worst case) 

la 

P 

0,0 

i 

P 2 

b 

T( s) 
a 0 

0.0828 

14.02 

86.54 

0.00686 

0.845 

2805 

2.01x 10~ 4 

3xl0~ 4 

0.19 

0.07 

4x 10~ 5 

0.005 

10 

O.lOx 10~ 4 

6xl0 -4 

0.44 

0.15 

l.OxlO" 4 

0.011 

24 

8x 10 -4 

0.59 

0.21 

1.4x 10 -4 

0.015 

33 


P 

0.0831 

2xl0~ 4 

X 

o 

1 

4- 

6x 10 -4 


Og 

13.70 

0.16 

0.32 

0.41 


i 

86.47 

0.06 

0.12 

0.15 

lb 

p 2 

0.00690 

4xl0~ 5 

OO 

X 

o 

1 

ai 

l.Ox 10 -4 


b 

0.844 

0.004 

0.008 

0.010 


T(s) 

2884 

10 

20 

25 


°o 

1.92x 10~ 4 

O.lOx 10 -4 




P 

0.0830 

3xl0' 4 

Cn 

X 

o 

1 

5x 10 -4 


Oo 

13.99 

0.24 

0.35 

0.37 


i 

86.55 

0.08 

0.12 

0.13 

2a 

p 2 

0.00688 

5xl0' 5 

OO 

X 

o 

1 

or 

8x 10 -5 


b 

0.841 

0.006 

0.009 

0.010 


T(s) 

2850 

13 

19 

20 


°0 

2.58x 10 -4 

0.13x 10 -4 




P 

0.0836 

4xl0' 4 

6xl0 -4 

9x 10“ 4 


Oo 

13.24 

0.22 

0.37 

0.54 


i 

86.27 

0.08 

0.14 

0.20 

2b 

P 2 

0.00699 

6x 10~ 5 

l.OxlO" 4 

1.4x 10 -4 


b 

0.861 

0.005 

0.009 

0.013 


T(s) 

2832 

15 

26 

37 


Oo 

2.88x 10~ 4 

0.15x 10 -4 
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Table 6: Transit parameter estimates for the four observations, by fitting p as a free parameter, cio and i 
identical for the observations at the same wavelength, with ATLAS quadratic limb darkening coefficients 
(Tab. [2J . We report the partial error bars obtained by the residuals, the final error bars, and the worst case 


error bars (Eq. Q] Sec. 12.21 and App. IA.3I). 


Obs. number 

Parameters 

Best values 

l-cr errors 

(residual scatter only) 

l-cr errors 
(ICA) 

l-cr errors 
(ICA worst case) 


Pla 

0.0831 

2.5x 10 -4 

5xl0“ 4 

7x10 -4 ' 


Pib 

0.0840 

2.5x 10 -4 

5xl0“ 4 

7x10 -4 


do 

14.04 

0.13 

0.28 

0.37 


i 

86.58 

0.05 

0.10 

0.13 

la + lb 

Pla 

0.00690 

4x 10 -5 

9x10 -5 

l.lxlO” 4 


pit 

0.00706 

4x 10 -5 

9xl0“ 5 

1.2x 10 -4 


b 

0.836 

0.004 

0.007 

0.010 


T(s) 

2876 

8 

16 

21 


c 0 

2.05x 10 -4 

0.07x 10 -4 




P2a 

0.0838 

3x 10 -4 

5xl0“ 4 

6x10 -4 


P2b 

0.0838 

3x 10 -4 

5xl0“ 4 

6x10 -4 


do 

13.74 

0.17 

0.26 

0.33 


i 

86.47 

0.06 

0.10 

0.12 

2a + 2b 

Pla 

0.00702 

5x 10 -5 

8xl0“ 5 

l.lxlO -4 


P\b 

0.00702 

5x 10 -5 

8xl0“ 5 

l.lxlO” 4 


b 

0.845 

0.004 

0.007 

0.009 


T(s) 

2868 

10 

16 

20 


Co 

2.73x 10 -4 

O.lOx 10 -4 
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Table 7: Transit parameter estimates for the four observations, by fitting p as a free parameter, cio and i 
identical for the observations at the same wavelength, with PHOENIX quadratic limb darkening coefficients 
(Tab. [2]). We report the partial error bars obtained by the residuals, the final error bars, and the worst case 


error bars (Eq. Q] Sec. 12.21 and App. IA.3I). 


Obs. number 

Parameters 

Best values 

l-cr errors 

(residual scatter only) 

l-cr errors 
(ICA) 

l-cr errors 
(ICA worst case) 


Pla 

0.0825 

2x 10 -4 

5xl0“ 4 

6x10 -4 ' 


Plb 

0.0834 

2xl0“ 4 

5xl0“ 4 

6xl0“ 4 


do 

13.82 

0.12 

0.26 

0.34 


i 

86.49 

0.04 

0.09 

0.12 

la + lb 

Via 

0.00681 

4x 10 -5 

8x10 -5 

1.0xl0“ 4 


Plb 

0.00696 

4x 10 -5 

8xl0“ 5 

l.Ox 10 -4 


b 

0.845 

0.003 

0.007 

0.008 


T(s) 

2845 

7 

16 

21 


Co 

2.04x 10 -4 

0.07x 10 -4 




P2a 

0.0833 

3x10“ 4 

5xl0“ 4 

6x10 -4 


P2b 

0.0833 

3x 10“ 4 

5xl0“ 4 

6x10 -4 


do 

13.57 

0.16 

0.25 

0.31 


i 

86.40 

0.06 

0.09 

0.11 

2a + 2b 

P2a 

0.00694 

5x 10“ 5 

8xl0“ 5 

l.Ox 10 -4 


Plb 

0.00694 

5x 10“ 5 

8xl0“ 5 

l.Ox 10~ 4 


b 

0.852 

0.004 

0.006 

0.008 


T(s) 

2840 

10 

16 

20 


Co 

2.74x 10 -4 

O.lOx 10 -4 
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